Preprint typeset in .THEP style - HYPER VERSION IFT-UAM/CSIC-08-02, hep-ph/yymmnnn 



Primordial magnetic fields from preheating 
at the electroweak scale 



Andres Dfaz-Gil Juan Garcfa-Bellido Margarita Garcia Perez " and 
Antonio Gonzalez-Arroyo 

00 " Instituto de Fisica Teorica UAM/CSIC and Departamento de Fisica Teorica 

Universidad Autonoma de Madrid, E-28049-Madrid, Spain 
CSj Kavli Institute for Theoretical Physics, University of California Santa Barbara, 

>^ CA 93106-4030, 

£'-mail; andres.diazgil@uam.es, juan.garciabellido@uam.es, 
^ margarita.garcia@uam.es, antonio.gonzalez-arroyo@uam.es 

Abstract: We analyze the generation of helical magnetic fields during preheating 
i-C in a model of low-scale electroweak (EW) hybrid inflation. We show how the inhomo- 

I geneities in the Higgs field, resulting from tachyonic preheating after infiation, seed 

§^ the magnetic fields in a way analogous to that predicted by Vachaspati and Cornwall 

r^, in the context of the EW symmetry breaking. At this stage, the helical nature of the 

generated magnetic fields is linked to the non-trivial winding of the Higgs-field. We 
>- analyze non-perturbatively the evolution of these helical seeds through the highly 

non-linear stages of symmetry breaking (SB) and beyond. Electroweak SB occurs 
via the nucleation and growth of Higgs bubbles which squeeze the magnetic fields 
into string-like structures. The 1^-boson charge density clusters in lumps around the 
magnetic strings. After symmetry breaking, a detailed analysis of the magnetic field 
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Fourier spectrum shows two well differentiated components: a UV radiation tail at 
a temperature T ~ 0.23 mjj, slowly growing with time, and an IR peak associated 
to the helical magnetic fields, which seems to follow inverse cascade. The system 
enters a regime in which we observe that both the amplitude {pb/pew ~ 10~^) and 
the correlation length of the magnetic field grow linearly with time. During this 
stage of evolution we also observe a power-law growth in the helical susceptibility. 
These properties support the possibility that our scenario could provide the seeds 
eventually evolving into the microgauss fields observed today in galaxies and clusters 
of galaxies. 
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1. Introduction 

The origin of magnetic fields (MF) is one of the remaining mysteries in relativistic 
astrophysics and cosmology (for reviews see the list of references [1]-[18]). Magnetic 
fields play an important role in the evolution of the primordial plasma in the early 
universe (possibly also in cosmic phase transitions), in the propagation of cosmic 
rays in our galaxy, as well as in clusters of galaxies. They may influence galaxy 
formation and large scale structures, and they may generate a stochastic background 
of gravitational waves. The connection between magnetic fields and gravitational 
waves is particularly intriguing. Since MF induce an anisotropic stress tensor, this 
can act as a source of gravitational waves (see [19]). Large amplitude magnetic 
fields from primordial turbulence could induce a significant stochastic background of 
gravitational waves which could be seen by LIGO or BBO, with a specific spectral 
signature. 

Magnetic fields have been found on the scale of galaxies and clusters of galaxies 
with a magnitude of order the microgauss. There is even some evidence of their 
existence on the scale of superclusters (for a review on observational results see [20]). 
Summarizing the measured MF values on all scales L : 

• galaxies: 5 ^ 50 jiG at L < 1 kpc; i? ~ 5 — 10 /iG at L ~ 10 kpc. 

• clusters: i? ~ 1 yuG at L ~ 1 Mpc. 

• superclusters: B < 10"^ - 10"^ /iG at L ~ 1 - 50 Mpc. 

• CMB: B < 10-3 - 10-5 /iG at L > 100 Mpc. 

• Primordial nucleosynthesis: B < 10^^ G at T = 10^ K. 

where the last bound (BEN) comes from the modification that such a background 
would imply for the expansion rate of the universe at primordial nucleosynthesis, 
which would change the observed Helium abundance. 

The main difficulty in understanding the origin of magnetic fields is not in their 
amplitude (i.e. magnitude) but in its correlation scale, from galaxies to clusters to 
superclusters. The microgauss order of magnitude of present galactic MF could be 
explained easily from an amplification via a dynamo mechanism initiated by a tiny 
seed, with B ~ 10"^^ — 10"^° G (when taking into account gravitational collapse in a 
flat ACDM model). The explanation of the scale of the magnetic seed in this case is 
rather straightforward. The dynamo mechanism is an exponential mechanism which 
makes the MF amplitude increase a factor e at every turn of the object (typically 
a galaxy) with free charge and thus large electrical conductivity. Since the typical 
galaxy has made around 30 turns in their lifetime, the growth factor is e^^ = 10^^. 
Since we observe microgauss, we just need a seed i?sccd ~ 10"^^ G over a scale of 
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30 kpc. This is the MF after gravitational collapse. Typically a galaxy forms by 
gravitational collapse of a lump of matter the size of about a Mpc with density of 
order the critical density, and ends collapsing to a size of order 30 kpc and density 
Pgai ~ lO^pc- By flux conservation, the gravitational collapse amplification gives an 
extra factor 

(Pgal/Pc)'/' ~ 10\ 

which gives a seed -Bseed ~ 10^^^ G over a scale of 1 Mpc. This calculation was done 
assuming matter domination. If we consider a ACDM universe, then gravitational 
collapse amplification is greater and the seed can start with -Bseed ~ 10^'^'^ G over a 
scale of 1 Mpc. This is the minimal value required for a typical galaxy. 

The microgauss amplitude at cluster scales is more difficult to explain via a 
dynamo mechanism because it did not have as much time since its formation to 
build up from such a tiny seed, and the order of fractions of microgauss amplitude 
at supercluster scales is simply impossible to explain by dynamo mechanisms or 
gravitational collapse. In any case, even in the presence of dynamo amplification, an 
initial magnetic seed is required which is not provided by the dynamo mechanism 
itself. Theoretical models trying to account for the origin of the primordial seeds can 
be classified in two groups: 

• Astrophysical: Biermann battery in intergalactic shocks, stellar magnetic winds 
(like in our Sun), supernova explosions, galactic outfiows in the inter-galactic 
medium (IGM), quasar outfiows of magnetized plasma into the intra-cluster 
medium (ICM), see Refs. [9, 12, 16], and a recently suggested proposal in 
conjunction with high energy cosmic rays [21]. 

• Cosmological: Early universe phase transitions [22]- [33], magnetic helicity to- 
gether with the baryon asymmetry of the universe (BAU) at the electroweak 
(EW) transition [34]-[41], via hypercharge and hypermagnetic field generation 
before EW transition [42, 43], from second order cosmological perturbations 
from infiation [44]-[63], from preheating after infiation [64]-[67], etc. 

Moreover, MF have also been observed in quasars at redshift 2; ~ 2, again with a 
magnitude of order the microgauss. This indicates not only ubiquity but also invari- 
ance (within an order of magnitude) with time. Such features cry for a cosmological, 
rather than astrophysical, origin of MF. Could it be that some yet unknown mech- 
anism directly generated microgauss MF on all scales? The first reaction is to ask 
about the dynamo mechanism in galaxies, would it not amplify this microgauss MF 
to even larger amplitudes, as can be seen in neutron stars, and even our Sun? The 
surprising answer is no, because a few microgauss is the maximum magnetic field 
possible on galactic scales, due to the existence of relativistic cosmic rays and ion- 
ized gas moving at large speeds. If one computes the total energy density in cosmic 
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rays (integrating the measured flux spectrum over aU energies), one finds 

\^PgWI(? = 0.5 eV/cm^ , 

and a similar number for the energy density in the ionized gas moving with rotation 
speeds of order 200 km/s, 

^Pgas'y^/c^ = 0.3 eV/cm^ . 

If we assume that magnetic fields are in equilibrium, due to their interaction with 
the cosmic rays and the gas, and furthermore we suppose equipartition, then their 
energy density (using 1 G = 1.95 x 10~^° GeV^) becomes 

Pb = B^/{8ti) = 0.5 eV/cia = (5 /xG)7(87r) , 

which corresponds to a few microgauss, in surprising agreement with observations. 
Some people suggest that this argument may also explain the cluster MF value. 

The ubiquity of MF with similar amplitude on all scales reminds us of the issue 
of Helium abundance in the universe. Early measurements in the fourties indicated 
that the Helium mass fraction to Hydrogen in the Universe was about a quarter, 
very nearly everywhere. This observation was correctly interpreted by Gamow and 
collaborators as indicating a primordial origin. Simple order of magnitude compu- 
tation of nuclear interaction rates (mainly those of deuterium, a necessary step in 
the reactions from H to He) and comparison with the rate of expansion in the early 
universe at temperatures of order the nuclear transitions (i.e. MeV), together with 
the then largely unknown neutron decay rate, suggested that the present abundance 
of Helium could have been produced from Hydrogen in the early universe and thus 
be present everywhere. The other light elements seemed to require further synthesis 
in stars and thus depended on location, but the Helium was ubiquitous because it 
was there from the very beginning. 

Something similar may have happened with magnetic fields, if they were gener- 
ated in the early universe by some unknown mechanism and then redshifted until 
today. The question is what is the typical energy density which today gives the 
order microgauss fields? These fields (if homogeneous) redshift as radiation, i.e. 
Pb{(i) = PB(today)(ao/a)^. Like with Helium, we have to ask what was the energy 
scale of interactions responsible for the generation of primordial magnetic fields? 
Photons are massless so in principle any scale, as long as there are charged particles, 
is sufficient to generate magnetic fields, and this is the reason why there is still so 
much debate as to their origin. However, was the universe always permeated with 
electromagnetic waves? The answer is no, the electromagnetic interaction as we 
know it came into being at a very precise time, when the electroweak (EW) force 
broke into the weak interactions plus electromagnetism. Before we could not talk 
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about photons and magnetic fields. This occurred when the typical energy (or tem- 
perature) in the universe was around Tew ~ 100 GeV. If we construct an energy 
density with this scale we get pew ~ 10® GeV*. At that time the universe was (or 
became) radiation dominated. If we now redshift this MF energy density until today 
(To = 2.725 K) we get 

PB(today) = (To/Tew)Vew ~ 3.04 x 10"^^ GeV* = 0.4 eV/cia 

which is precisely the order of magnitude of the present MF energy density.^ This 
would be enough to explain the cluster and supercluster values, and would perhaps 
require a mild dynamo mechanism to grow to galactic values (if the fraction / ^ 1). 
The question is whether this is just a coincidence or it is hinting directly at its 
origin.^ While other mechanisms require a seed with an arbitrary scale (typically 
B ~ 10~^^ G, so that today we observe microgauss MF on galactic scales via the 
dynamo mechanism), there is no physical reason behind this scale. On the other 
hand, the EW scale is a natural scale for the generation on magnetic fields since 
it is the scale at which electromagnetism arises for the first time as a fundamental 
interaction. 

Whether this is sufficient reason to assign the EW energy scale to the origin 
of magnetic fields is another issue. In particular, it is not clear how to obtain the 
large correlation length of magnetic fields observed at galactic and cluster scales. 
Any physical mechanism that creates magnetic fields must be necessarily causal, but 
at high temperatures in the early universe there is also a natural coherence scale 
given by the particle horizon. At the electroweak scale the physical horizon is 10~^° 
light-seconds (~ 3 cm), which today corresponds to a co-moving scale of 0.3 mpc 
(~ 1 AU), clearly insufficient when compared even with the irregular (turbulent) 
component of the galactic magnetic field {L ~ 100 pc), not to mention the regular 
(uniform) component, which has correlations L ~ 10 kpc. It thus seems impossible 
to explain the coherent magnetic fields observed on galaxy clusters and supercluster 
scales (of order 10 Mpc) with intensities of order fiG to nG. 

There is however a second coincidence, which makes things even more intriguing. 
If we assume that the plasma after the electroweak transition is sufficiently turbulent 
to maintain magnetic fields of the largest possible coherence scales via inverse cascade 
[29]-[33], then we could reach cosmological scales today. Let us follow the argument. 
The largest coherence scale at the electroweak transition is the physical horizon, of 
order 3cm. If a strong inverse cascade is active, then the coherence length of the 
magnetic fields will grow as fast as the horizon (it cannot grow faster). This means 

""^We could be even more conservative and suppose that the fraction of magnetic field energy 
density to radiation at the time of the EW transition was given by / = Ps/Piad < 1- In this case, 
the present MF magnitude would be Bq ~ 5 f^^^ iiG. 

^Some authors suppose that the generation occurred earlier in the form of hypcrmagnctic fields 
and was then converted into ordinary magnetic fields at the EW scale [43]. 
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that it grows like the scale factor squared during the radiation dominated era. This 
ideal situation could only last while there is a plasma and thus it is bound to stop 
acting at photon decoupling, when the universe becomes neutral. Since then, the 
correlation length can only grow with the expansion of the universe, as the scale 
factor. If we take this effect into account from the electroweak scale until today we 
find, using the adiabatic expansion relation T oc a~^, 

eo = ^Ew (—) ' ^ = 3 cm ' ^ ~ 6 X 10^^ cm = 20 Mpc , (1.1) 

where we have made the approximation that equality and decoupling occurred more 
or less simultaneously (a careful computation gives only a minor correction). The 
surprising thing is that this simple calculation gives precisely the order of magnitude 
for the largest correlation length of cosmic MF ever observed (i.e. cluster scales). If 
the agreement in the magnitude of the primordial MF seed seemed peculiar, the fact 
that an inverse cascade could also be responsible for the observed correlation length 
becomes a surprising coincidence, probably hinting at an underlying mechanism. 
It is therefore worthwhile exploring the conditions that could have taken place at 
the electroweak transition which could give rise to a significant fraction of energy 
density in magnetic fields, and be responsible for a sustained period of inverse cascade 
until photon decoupling. It has been shown in Refs. [30]- [33] that one important 
ingredient is the generation of magnetic fields with a non trivial helical component, 
which guarantees an optimal amplification of the magnetic correlation length through 
inverse cascade. A very good account of the large number of works investigating these 
issues, with a complete list of references is given in Ref. [9] (see also [27]- [43]). 

In this paper we propose a scenario in which the electroweak transition takes 
place at the end of a brief period of hybrid inflation^ It has been conjectured that 
preheating and early reheating in this model could provide an alternative mechanism 
to generate the baryon asymmetry in the universe [70]- [76] and a way to source 
gravitational waves [77]. In this paper we analyze whether it could also give rise to 
primordial magnetic fields with the required amplitude and correlation length. This 
issue has been partially addressed in a recent letter [78]. Here we will present a 
complete account of the results obtained and a detailed description of the approach 
employed in the analysis. Our set up provides a specific realization of some of the 
proposals described above. In particular, we will see how helical magnetic fields arise 
from the inhomogeneities in the spatial distribution of the Higgs field, along the lines 
conjectured by Vachaspati [34, 36] and Cornwall [38] some years ago. 

The paper is organized as follows. In section 2 we describe the hybrid inflation 
model that we will be using and revise, following Ref. [73], how to solve the quan- 

^Note that we do not need a 60 e-fold period of inflation, just a few (~ 5) e-folds of low 
scale thermal inflation [68, 69] to cool down the universe. The amplitude of CMB temperature 
fluctuations would be determined by the usual 55 c-folds of high-scale (e.g. GUT) inflation. 
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turn evolution of the system from the end of inflation until non-linearities start to 
become important. Beyond this time, a fully non-perturbative approach is required. 
Fortunately, the time evolution at this stage can be described within the classical ap- 
proximation as demonstrated in Ref. [73]. Details on the methodological set up and 
the lattice implementation are presented in section 3. Section 4 analyzes the mech- 
anism leading to helical magnetic fields following from the inhomogeneities in the 
Higgs field, which are seeded by the Higgs quantum fluctuations that arise from the 
period of linear quantum evolution. Strings of magnetic flux, carrying non- vanishing 
helicity, are clearly observed. They persist and are even enhanced as the system pro- 
gresses towards the true vacuum. We also analyze here the structure of the plasma 
of H^— charges which accompany the magnetic fields during this period. The fate 
of these magnetic fields at later times is discussed in section 5, where we present a 
detailed study of the spectrum of the magnetic field. We will show that there is a 
significant helical magnetic field remnant whose amplitude and correlation length are 
amplified linearly in time. In section 6 we will discuss lattice and finite volume inde- 
pendence of our results, as well as the dependence of magnetic field production on the 
Higgs- to W-mass ratio. Conclusions and prospects for further work are presented 
in section 7. Finally, a few technical points about the lattice discretization of the 
classical equations of motion, the Maxwell equations and electromagnetic radiation 
are described in the Appendices A, B and C. Appendix D is devoted to an analysis 
of the Gaussian random fields that provide the initial Higgs field distribution. 

2. The model 

The scenario we will be considering is that of preheating after a period of hybrid 
inflation which ends at the EW scale. This was first introduced in Ref. [70] to 
provide a new mechanism for the generation of baryon asymmetry in the Universe 
(BAU). It has been extensively studied since then both in connection with BAU [70]- 
[76] and in relation with the production of gravitational waves [77]. In this paper 
and in Ref. [78] we include for the first time the Hypercharge field in order to study 
the generation of electromagnetic fields during preheating (preliminary results can 
be found in [79]). In this section we will introduce the model and describe the first 
stages of evolution after inflation ends which provide the initial conditions for the 
non-linear approach addressed in section 3.1. 

The Hybrid inflation model is attained by extending the Standard Model with 
the addition of a scalar field, the inflaton, singlet under the gauge group. The scalar 
sector thus includes the Higgs field: $ = ^(0° H + i(f)"'Ta) {ja are the Pauli matrices) 
and the singlet inflaton x which couples only to the Higgs via the scalar potential: 

V($, x) = Vo + \{g\' - m^) \<\>f + ^10^ + \^\^ , (2.1) 
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where |0|2 = 2Tr<l>t$, jjL is the inflaton mass in the false vacuum and = \/2 m = 
\/2Xv the Higgs mass, with ^ = 246 GeV the Higgs vacuum expectation value at the 
electroweak scale. The gauge sector contains both the SU(2) and the hypercharge 
U(l) fields with 



abc 



and 



■Y 



d B 

their respective field strengths. The covariant derivative is: 



% 1 



(2.2) 
(2.3) 

(2.4) 



with g^ the SU(2) gauge coupling and g^ the hypercharge coupling. In this work we 
can safely ignore fermionic fields since the time scales involved in the perturbative 
decay of the Higgs field into fermions are much larger than the ones considered here. 
With all these definitions the Lagrangian density of the model becomes: 



y($,x). (2.5) 



For our analysis we have fixed the W mass and the Z to W mass ratio to the 
experimental values [80]. We have analyzed three different values of the Higgs to W 
mass ratio: m^/m^ = 2\/2\/ g^ = 2, 3 and 4.65. The Higgs-infiaton coupling has 
been fixed to g^ = 2A as in super-symmetric models [71, 81] and we have taken the 
infiaton bare mass fi = 10~^gv ~ 0. 

The extraction of the electromagnetic content of the SU(2)xU(l) fields in the 
Lagrangian proceeds in the usual way. Fixing the unitary gauge, $(x) = p{x) U, the 
Z-boson field and the electromagnetic field are extracted from appropriate orthogonal 
combinations of the SU(2) and hypercharge vector potentials: 



cos 9w AJx) 



sin 6w BJx) 



(2.6) 
(2.7) 



with <f{x) = $(x)(l,0)^ the Higgs doublet. This separation can only be done unam- 
biguously when the Higgs field is on the true vacuum, i.e. in the broken symmetry 
phase. However even in that phase there can be points where the Higgs field vanishes 
and the symmetry is locally restored (a typical example of a configuration exhibiting 
such behavior is the sphaleron). At those points there is no unique way to define the 
electromagnetic fields. In fact 't Hooft was the first to point out in Ref. [82] the con- 
sequences of this ambiguity in the Georgi-Glashow model, tying it to the appearance 
of non-trivial configurations like monopoles or strings, acting as sources of magnetic 
fields. In Ref. [34] Vachaspati pointed out that a similar mechanism was at work in 
the electroweak model where the sources for magnetic field generation are tied to the 
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presence of non-homogeneous phases in the Higgs field. In the following sections we 
will analyze in detail how this mechanism is realized during the period of preheating 
after infiation. 

2.1 Linear quantum evolution 

Following Refs. [73, 74], we will address here the first stages of evolution starting at 
the end of infiation. 

The period of infiation is characterized by the fact that the Higgs and infiaton 
fields are displaced from the true minimum of the potential. In this case, infiation 
is driven by the false vacuum energy, Vq = During this time the infiaton 

homogeneous mode, xo = {x)j dominates the dynamics. After only about 5-10 e- 
folds the Universe has cooled down and all other particle species have been diluted, 
remaining in the de Sitter vacuum.^ 

The interaction between the Higgs and infiaton fields drives the end of infiation 
and triggers EW symmetry breaking. The way this proceeds is as follows. Close to 
the time when infiation ends, denoted by tc, the time evolution of the infiaton zero 
mode can be approximated by: 



where Xc = Xo{tc) = i^/g- Here V denotes the infiaton dimensionless velocity, 
defined through this equation and fixed to = 0.024 in our analysis [73]. The 
variation of Xo(^) induces, via the Higgs- infiaton coupling, a time dependence of the 
effective Higgs mass parameter, = —M^{t — tc) = —2VnT'{t — tc)-, which changes 
from positive to negative, triggering electroweak symmetry breaking. Accordingly, 
the time when infiation ends, tc? is characterized as the critical point where the Higgs 
field becomes massless. 

As described in detail in Refs. [73, 74], it is possible to solve exactly the quantum 
evolution of the system around tc if non-linearities in the Higgs field and the inter- 
action with the gauge fields are neglected. As we will see below, this is a reasonable 
assumption at this stage. In this approximation the Higgs field is effectively described 
as a free scalar field with a time dependent mass m^{t). Its quantum evolution can be 
solved in terms of Airy functions [73]. After tc, low momentum modes of the Higgs 
field grow exponentially in a process known as "tachyonic preheating" [71]. Due 
to the tachyonic growth, low momentum Higgs field modes acquire large occupation 
numbers and, soon after tc, they evolve as classical modes. This is a very fast process 
so that all other modes can be taken to remain in the quantum vacuum (ground) 
state, justifying thus the linear approximation. These modes will be later populated 
through the interaction with the Higgs field tachyonic modes. Once non-linearities 

^For clcctroweak-scale inflation and the range of momenta we will be considering, de Sitter 
vacuum is indistinguishable from the Minkowski vacuum. 
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start to become relevant the approximation ceases to be valid and a full non-linear 
treatment is required. Our strategy for dealing with the later stages of evolution will 
be presented in the next section. 



3. Methodological set up 

Beyond the quantum linear evolution described in the previous section we have to 
deal with the non-linear dynamics of the Higgs field and its coupling to the gauge 
fields. Our approach is based upon the classical approximation (details can be found 
in [73, 74] - see also [76]). The validity of this approximation relies on the fast growth 
of tachyonic modes as explained previously. In what follows we will describe several 
aspects of our procedure. 

3.1 Initial conditions for the the non-linear evolution 

As mentioned previously the initial stages after the end of inflation (t = tc) lead 
to a rapid growth of the tachyonic modes which tend to behave classically. The 
correlation functions of the Higgs field resulting from the initial quantum evolution 
can be computed. Our approach is to use these results as initial conditions for the 
classical evolution of the system. The quantum fluctuations translate into stochastic 
initial conditions for the Higgs field, whose correlations are designed to match the 
Weyl-ordered quantum expectation values. The matching of the two methods is done 
at an initial time t = ti > tc that must be large enough for classical behaviour to set 
in and small enough to make the non-linear terms small. This leaves a window of 
possible values of t,. We tested the robustness of the results with respect to changes 
in ti within these limits, giving confidence on the self-consistency of our approach. 

Given the linear character of the initial quantum evolution, the Higgs field mo- 
mentum modes (j)'^, at t = ti behave as Gaussian random variables of zero- mean 
following a Rayleigh distribution: 



where 6'^ is the phase of the complex random variable 0^. The dispersion of the 
modulus is expressed in terms of the power spectrum P{k,ti) = P((T^)^, and can 
be computed analytically in terms of Airy functions [73]. For practical purposes it 
is better to work with a simple functional fit to the power spectrum (Eq. (D.l)). 
Notice that we have introduced a momentum cut-off, removing modes which have 
not become tachyonic (classical). As explained in Ref. [73], this is compensated by 
a renormalization of the parameters. 

The study of the properties of this 4-component Gaussian random field is col- 
lected in Appendix D. Its features depend on several parameters: the Higgs mass, 
the initial infiaton velocity V, the momentum cut-off and the choice of initial time 




(3.1) 
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tj. The first two appear combined in a new scale M = (2VY^^m characteristic of the 
initial linear evolution. The main conclusion drawn in Appendix D, is that to a large 
extent all these parameter dependencies translate into setting two main scales: a spa- 
tial length scale, ^o, and the Higgs dispersion, a, which determines the magnitude of 
the Higgs field. 

More specifically, in Appendix D we study the distribution of local maxima in 
|0(x)|. These "peaks" are the seeds that will later grow with time and develop into 
bubbles which start expanding and colliding among themselves once the Higgs fields 
enters the non-linear regime characteristic of symmetry breaking. This process was 
described in detail in Ref . [73] . Note that the multicomponent character of the Higgs 
field affects the results but, more importantly, it gives rise to new observables, some 
of which are intimately connected to the physical phenomena which are the main 
goal of this paper. This will be described in the next section. 

To complete the description of the initial conditions, we mention that, similarly 
to the high-momentum modes of the Higgs field, all other non-tachyonic modes are set 
to zero. These include the non-homogeneous modes of the inflaton and the vector 
potentials of the SU(2) and hypercharge gauge fields. The initial time-derivatives 
of these quantities are also set to zero except for the gauge fields which have to 
be chosen such that the Gauss constraint is satisfied as an initial condition. The 
dynamic equations guarantee that the constraint will continue to hold at later times. 
The aforementioned robustness of the results to the choice of initial time ti implies 
that our physical conclusions do not depend on minor modifications of these initial 
conditions. 

3.2 Numerical procedure 

In order to study the non-linear evolution of the system with our stochastic ini- 
tial conditions we have made use of the lattice approach. This has the advantage 
that classical equations of motion are discretized preserving full gauge invariance of 
the system. Generally speaking the procedure is standard. Details on the lattice 
Lagrangian and the lattice form of the equations of motion are presented in Ap- 
pendix A. At early times the errors associated to discretization are very small due 
to the cut-off form of the initial spectrum. This shows up in the very mild depen- 
dence of the results on the spatial a and temporal at lattice spacings. This contrasts 
with other situations in which lattice techniques have been used. As time evolves 
higher momenta of the fields grow and start to play a role, eventually leading to 
a breakdown of the approximation. We have explicitly analysed that this does not 
occur for the range of times covered in this paper. A different approach needs to be 
followed if one wishes to reach times in which full thermalisation has been reached. 
Notice, however, that this goal also demands the introduction of fermionic degrees 
of freedom which can be safely ignored in our time span. Our present results can be 
used as initial conditions for the study of the late time behaviour of the system. 
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N 


777 n 








65 


1/40 


150 


80 


0.52 


1/40 


0.150 


100 


0.42 


1/40 


0.150 


100 


0.52 


1/40 


0.125 


100 


0.65 


1/40 


0.100 


120 


0.65 


1/40 


0.080 



Another approximation needed 
for the numerical procedure is to 
put the system in a box with pe- 
riodic boundary conditions. The 
physical volume, V = L'^ is given 
in terms of the minimum momen- 
tum: L = 27r/p^i„. The latter has 
to be chosen judiciously to lie well 
within the tachyonic band of the 
Higgs field. The dependencies of 
the results can be monitored by 
using different values for the pa- 
rameters of the simulation. In Ta- 
ble 1 we enumerate the different 
lattice sizes, spacings and physi- 
cal volumes that we have used. 

Due to its relevance for the 
goals of our paper, we will now 

explain in detail how the electromagnetic and Z fields are defined in our lattice 
approach. This can only be done unambiguously when the Higgs field is in the true 
vacuum, i.e. in the broken symmetry phase. One can compute, in a gauge invariant 
way, the field associated to the Z-boson potential as: 



Table 1: List of lattice parameters: a and at 
are respectively the spatial and temporal lattice 
spacings, Ng is the number of lattice points and 
Pmin = 271 /{NgO) is the minimum momentum. The 
Ns = 120 lattice has only been used for the study 
of the initial configuration. The number of different 
configurations of each lattice ranges from 80 to 200, 
depending of the lattice and the choice of parame- 
ters. 



-iTr 



n(D^<l>(m))$t( 



-iTr 



(m)||0(m + 



'|0(m)| 



m-h/i)| 



(3.2) 

ai,9zZ^{x) , (3.3) 



where we have introduced the adjoint unit vector h = UaTa, with components: 



na{x) 



(3.4) 



with (f{x) 
gz and 



<l>(x)(l,0)^ the Higgs doublet. The Z boson coupling is denoted by 
Ct, (Xq = (It. is the lattice covariant derivative operator defined 



in Eq. (A. 6) of Appendix A. 



and B^{n) 



are, respectively, the SU(2) and 



hypercharge link fields introduced in Appendix A. Notice that continuum quantities 
are defined with calligraphic letters to distinguish them from the lattice quantities. 
Our definition of the Z boson potential corresponds to the standard one in the unitary 
gauge. 
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We define the Z boson and hypercharge field strengths through the clover aver- 
ages ^ 

F^^{'m) = {Af,Z^{m)-A^Zf,{m))ciov ^ af,a^gzJ^^„{x) + 0{a'^) (3.5) 

and 

F^^{'m) = {Af,e„{m)-AJ^{m))c\ov ^ af,a,yg^T]^^{x) + Oia"^) , (3.6) 

where B^{m) = exp(i6'^(m)r3/2) is the hypercharge link, is the lattice derivative 
operator introduced in Eq. (A. 8) and (O)ciov denotes the clover averages defined in 
Eqs. (B.2)-(B.3). In terms of them we can compute the lattice electromagnetic field 
strength as: 

F;,(n) = sm\ew)F^,{n) - Fj,{n) a,a,eJ^;,{x) + 0{a') , (3.7) 

where jFJj, is the corresponding continuum electromagnetic field strength. This pro- 
vides a lattice gauge invariant definition of the electromagnetic field which is equiv- 
alent to the usual definition in the unitary gauge. 

4. The mechanism underlying magnetic field generation 

In this section we study the production of magnetic fields during the first stages of 
our EW preheating scenario. This analysis is performed in two steps. The first is to 
investigate the presence, size and structure of the magnetic fields generated by our 
Gaussian random field initial distribution. This complements the results presented 
in Appendix D. Then we will track the evolution of these magnetic fields through 
the highly non-linear stages associated to EW symmetry breaking. This is a crucial 
period where there are no viable alternatives to our methodological approach. 

4.1 Initial Magnetic fields 

A close look at our expression of the photon field reveals that Abelian electromag- 
netic fields are present in the first stages of the evolution. The discussion on how this 
comes about follows a line of argument very similar to that developed by Vachaspati 
in Ref. [34]. The tachyonic preheating phase leads to a multicomponent Gaussian 
Higgs field. The SU(2) and hypercharge gauge fields remain very small. This is incor- 
porated into our initial conditions by setting the hypercharge and SU(2) magnetic- 
like fields to zero and fixing the corresponding electric fields in order to satisfy the 
Gauss constraint. We work in the = gauge, which on the lattice corresponds 
to Un{t = ti) = Bn{t = ti) = 11. Projecting onto the Z and electromagnetic fields we 

^Corrections to the continuum approach of the time-space clover averages are order O{aoa^). 




Figure 1: (Left) Locus of points where the value of the Higgs field norm is below 0.03 f . 
(Right) Locus of points where twice the magnetic energy density (|-Bp) is above 0.03 m^. 
Data correspond to mt = 5 and = 2m^. 



obtain: 



Tr 



Tr 



9z 



h^l{x)d^il\x) 

n I dyQ{x)d^ri\x) — d^^l{x)dui^\x 

h [duVt{x)d^Vt\x 



(4.1) 



9w 



d^il{x)d^il\x] 



expressed in terms of the SU(2) matrix: 



Q(x) 



$(x) 



(4.2) 



It becomes clear that electromagnetic fields are sourced by the presence of inhomo- 
geneities in the Higgs field orientation. This is one of the essential ingredients in 
Vachaspati's proposal for magnetogenesis. 

The size and spatial distribution of this initial electromagnetic and Z fields can 
be obtained from the multicomponent Gaussian random field. In Appendix D we 
displayed the histogram of magnetic field values. Here we will focus on another aspect 
which is particularly interesting for the later evolution. This is the spatial distribution 
of points where the magnetic field intensity is larger. To investigate this, we show 
in Fig. 1 a 3-dimensional plot displaying the locus of points where the magnetic 
energy density is above 0.03 m^ for our initial configuration at mtj = 5. Notice 
that the regions of higher magnetic energy density exhibit a string-like geometry. 
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1200 




Figure 2: Histogram of Higgs winding number for the initial configuration mt = 5 for 
Pn,in = 0.15 m. 

Indeed, this spatial distribution tracks the location of regions of low Higgs field 
value, which are also presented in the figure. Although, the strings seem to end at 
certain spatial points, this is simply a reflection of the spreading of magnetic flux 
lines. Our electromagnetic field satisfies the Maxwell equations without magnetic 
sources or sinks. According to our formulas the initial magnetic-like component of 
the Z-boson field strength is directly proportional to the electromagnetic field and 
has identical structure. 

There is another important feature of magnetic fields which we have investigated. 
It corresponds to whether the initial field gives rise to a sizable helicity. In a finite 
volume the total magnetic helicity is defined both in configuration and in momentum 
space as: 



H 



d^x h(x) 



(fxA- B 



—I 



Y,Tfr-{my<B\k)), (4.3) 



where V is the volume of space. Notice that this equality makes use of Maxwell's 
condition V-B = 0, which is ensured by our magnetic field definition (3.7). At our 
initial time, by virtue of Eqs. (4.1), this quantity is proportional to the winding of the 
Higgs field. This is defined as the index of the map from the spatial volume to the 
group SU(2)=S'3, provided by the matrix VL{x). A histogram of the winding obtained 
for our initial Gaussian random field configurations is displayed in Fig. 2. The 
data are well described by a Gaussian distribution. Since we have not included CP 
violating terms, the mean value of the winding number is zero. However, we observe 
a non-zero dispersion from which one can obtain a non-zero volume-independent 



topological susceptibility x = 0.52 x 10" 



This translates into a corresponding 
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mt 



Figure 3: We show the time evolution of the electric (transverse Et and longitudinal El) 
and magnetic energy densities averaged over 150 configurations for m^^ = 3m^, ma = 0.42 
and Pmin = 0.15 m. 

non-vanishing helical mag netic susceptibility xh = (^^)/V = 0.38(3) m^. The Z 
helical susceptibility at this initial stage is xz = tan~^ 9w xh- 

In the next subsection we will study the evolution of this helical magnetic field 
during the highly non-linear epoch of symmetry breaking. This provides a connec- 
tion between magnetic field helicity, Z-strings and the occurrence of configurations 
carrying non-trivial Chern-Simons number. This result, which relates baryon num- 
ber generation and magnetic helicity, has been proposed, although along somewhat 
different lines, by Cornwall [38]. The connection has been studied recently by Copi et 
al. [83]. They showed that the sphaleron decay indeed gives rise to helical magnetic 
fields. 

At later stages, the Chern-Simons number creation processes stop, leaving behind 
a remnant magnetic helicity component [36], which is preserved in a plasma with 
high electrical conductivity. Thus, this could provide a signature of EW generation 
of primordial magnetic fields. 

4.2 Magnetic strings through symmetry breaking 

We will now focus upon the evolution of the system from the initial Gaussian random 
field situation until the onset of symmetry breaking. To have a global picture of the 
process we show in Fig. 3 the time evolution of the expectation value of the Higgs 
field from the initial time mti = 5 of our classical evolution. Notice the strong initial 
oscillations for times smaller than mt = 20, which are then progressively damped at 
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L L 

' X ' X 

Figure 4: Top: (Left) The location of tlie bubbles in the Higgs field norm (in red) with 
a lower cutoff set at 0.7 v and the locus of points with twice the magnetic energy density 
(|i?(x)P) (in blue) higher than 0.01 m^. (Right) Locus of points where the magnetic energy 
density is above 0.03 m^. Bottom: (Left) Two-dimensional contour plots of the Higgs 
field norm. (Right) Two-dimensional contour plots of the magnetic energy density. Data 
correspond to mt = 15 and = 2m^. 

larger times. The figure also displays the fraction of the total energy density carried 
by electromagnetic fields. We split it into its magnetic and electric components, 
and for the latter we analyze separately longitudinal and transverse parts. ^ We 
observe that between mt = 10 and mt = 15, there is an explosive growth of the 
electromagnetic fields correlated with the first minimum in the oscillation of the 

^The technicalities involved in the lattice definition of transverse and longitudinal fields as in 
the definition of the W bosons charge densities and currents are discussed in Appendix B. 
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1 z 

Figure 5: Top: (Left) Helicity of the magnetic field. (Right) Hehcity of the Z-boson field. 
Bottom: (Left) Two dimensional contom' plots of the helicity of the magnetic field. (Right) 
Two dimensional contour plots of the helicity of the Z-boson field. Data correspond to 
mt = 15, for m-jj = 2m^. 



Higgs field expectation value. The data in the figure corresponds to = 3m^, but 
similar behaviour is observed for the other ratios studied. 

We will now present the spatial structure observed for the magnetic fields at 
mt = 15 after the strong oscillation region. The corresponding distribution of the 
Higgs field modulus has been presented in Ref . [73] . There we showed that the initial 
Gaussian peaks lead to bubbles which expand and collide with neighbouring ones. 
This is illustrated in the top left of Fig. 4 where we display a snapshot of the Higgs 
field norm at mt = 15. At this time bubble shells (in red), that have grown out 
of the peaks in the initial Gaussian random field, fill almost all the volume of the 
box. Magnetic fields (shown in blue in the figure) appear as string-like structures 
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Figure 6: Averaged PHi{rQ), Eq. (4.4), at mt = 15. We also show the P asymptotic 
behaviour for the m^/m^ = 4.65 model. The other data corresponds to m^/m^ = 2. 

localized in the region between bubbles, where the Higgs field remains closer to the 
false vacuum for a longer period of time. This linkage between magnetic strings 
and Higgs field minima is even more evident in the two dimensional contour plots 
presented in the bottom half of Fig. 4. 

The structures observed in the regions of maximal magnetic density are repro- 
duced when looking at the helical part alone. This is exemplified by the comparison 
of Fig. 5 with Fig. 4. The figure also shows how the correlation between magnetic 
and Z boson fields, implicit in our initial conditions, is still preserved once gauge 
fields and non-linearities have started to play a role. An interesting observation can 
be made here concerning the connection with baryon number generation. Analysis 
of the cold EW transition show that sphaleron-like configurations, with non-trivial 
Chern-Simons number, are also located between bubble shells [74]- [75]. For non- 
zero Weinberg angle, sphalerons look like magnetic dipoles [84] and it is tempting to 
correlate the observed helical magnetic flux tubes with the alignment of sphaleron 
dipoles. Although a detailed investigation of this correlation is beyond the scope 
of this paper, our results for the distribution of magnetic helicity do indeed hint in 
that direction. An evaluation of the net helicity at late times and a discussion on its 
persistence will be postponed to section 5. 

In the previous figures, the closed string-like structure of the helicity and mag- 
netic field appears much more clearly that in the Gaussian random field initial condi- 
tion at mti = 5. To quantify the string- like character, we have analyzed the following 
quantity: 




(4.4) 



where h{x) denotes the hehcity density and the integration is on a box of length /, cen- 
tered at a point Tq at the center of one of the strings. Figure 6 shows the /-dependence 
of /^ifi(ro), averaged over several configurations. The figure is intended to show the 
one-dimensional character of the distribution in accordance to our string picture. In 
that case, Hiiro) should be /-independent in contrast with the /^-behaviour charac- 
teristic of an isotropic distribution. Both regimes are clearly observed in the figure. 
The stringy behaviour is displayed up to ml ~ 10, beyond which the plot shows how 
the data tends nicely to a straight line of slope equal to 2. This is to be expected 
once the box is big enough to contain several strings. This leads to an estimate of 
the string separation of m-^J, ~ 14, which is a significant fraction of the total length 
of the box. 

4.3 Charge lumps around magnetic field lines 

Up to now we have focused on the distribution of magnetic and Z-boson fields, but 
there is important additional information on the nature of the primordial plasma 
during these stages of preheating. Note that our initial conditions provide a source 
for charged VT-currents and a non-trivial charge density. It turns out that there is 
charge separation at the initial stage. Positive and negative charges are clustered 
into separate lumps which track the magnetic field lines. Figures 7 and 8 show this 
effect at mt = 15 and mt = 10, respectively. Note that there is a strong correlation 
between the magnetic field lines and the distribution of charges of opposite sign 
around them. The effect is seen particularly clear at early times, mt = 10, where 
the magnetic flux tubes are well defined, and there are fewer of them. The charge 
separation is consistent with the effect that would be produced by a combination 
of the drift currents induced by gradient and curvature effects from the magnetic 
field. The electric field is also strongly correlated with the location of the charge 
lumps, as expected. This charge separation might be responsible for the very slow 
screening observed for the longitudinal electric field, which will be discussed in the 
next section. 

The plasma generated during the first stages of evolution is, as we have shown, 
somewhat different from standard MHD plasmas (composed mainly of protons and 
electrons, together with photons). Here, long range string-like structures are observed 
in the electromagnetic fields, and opposite IV-charges cluster in large regions of 
space inducing non-trivial electric fields. It is expected that these charge lumps 
will eventually disintegrate when the VT-fields decay into light fermions (quarks and 
leptons), which travel at the speed of light and diffuse the charge, leading at late 
times to a standard MHD plasma. 

5. Late time evolution 

In order to claim a mechanism for cosmological magnetogenesis, the essential question 
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Figure 7: Top Left: Locus of points with magnetic energy density above 0.01 m^. 

Top Right: Locus of points with electric energy density above 0.01 m^. Bottom: 

The distribution of charge density, tracking the magnetic field lines. Pink and blue 
areas represent negative and positive charge densities respectively. Data correspond to 
mt = 15, for TTijj = 2m^. 

is whether the amplitude and correlation length of the generated fields are enough 
to seed the large scale magnetic fields observed today. In this section we will present 
evidence that a significant fraction of long range helical magnetic fields remains after 
EW symmetry breaking and is even amplified at later times, a period in which kinetic 
turbulence has been observed [85, 79]. As we will see below, our estimate for the 
amplitude of the magnetic field seed gives a fraction ~ 10^^ of the total energy 
density at the EW scale. This could be enough to seed the cluster and supercluster 
values without the need for a dynamo mechanism. 

More difficult is to address the issue of whether the magnetic field spectrum 
experiences inverse cascade, i.e. transference of energy from high to low momentum 
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Figure 8: Left: Locus of points with magnetic field density above 0.01 m^. Right: 

The 2-dimensional VF-charge distribution locahzed in lumps of opposite sign facing each 
other. Note that the location of the charge lumps is strongly correlated with the magnetic 
field flux tubes. These figures correspond to early times, mt = 10, for m,-^ = 2>m^. 

modes [29]- [33]. Inverse cascade is required to make the coherence length of the mag- 
netic field grow (almost) as fast as the horizon until the time of photon decoupling. 
Our approach does not allow to extrapolate the time evolution for sufficiently long 
times. Nevertheless, we will provide some evidence that inverse cascade might be at 
work. However, additional work is required to analyze if it can be sustained for a 
sufficiently long time. This might require a full magnetohydrodynamics treatment 
of the time evolution for which our set up will provide an initial condition. 

5.1 Magnetic helicity and electromagnetic energy densities 

We will first analyze in detail how electromagnetic fields evolve in time, paying 
particular attention to the evolution of the magnetic field helicity long after SSB. 

As mentioned above, the relevant quantity for helicity in the absence of CP vi- 
olation is the helical susceptibility xh- Its time evolution, for different values of 
the m^/m^ ratio, is displayed in Fig. 9. At the same time we display the helical 
susceptibility of the Z boson magnetic field, rescaled by tan'^ 9w to make it agree 
with the initial electromagnetic helicity, see discussion after Eq. ( 4.3). The late time 
behaviour, after mt ~ 60, gives further support to the Vachaspati-Cornwall's conjec- 
ture. It corroborates that, while the Z boson helicity is damped in time, the magnetic 
helicity is preserved and even increases with a power law dependence in time given 
by with a= 0.7(1), 0.8(1), 0.3(1) for m^/m^ = 2,3 and 4.65 respectively. The 
corresponding helical susceptibilities at mt = 100 are 0.11(2), 0.26(1), 0.12(2) m^. 
Note that the model with = 3m^ is more efficient than the others in generating 
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Figure 9: We display the time evolution of the helical susceptibilities for the magnetic 
field (with fit t") and the Z-boson field (with fit t^). The latter is rescaled by tan'^6\y to 
match the initial electromagnetic helicity. Top left is for m^lm^ = 2, averaged over 80 
configurations, with a = 0.7(1) and /5 = —0.27(4). Top right is for m^jm^ = 4.65, aver- 
aged over 80 configurations, with a = 0.3(1) and /? = —0.33(5). Bottom is for m^/m^ = 3, 
averaged over 200 configurations, with a = 0.8(1) and (3 = —0.82(4). All data correspond 
to ma = 0.42 and p^^^-, = 0.15 m. The top left figure also shows the time evolution of the 
Higgs mean to illustrate the time when SSB takes place. 



helicity at late times. This suggests a non monotonic dependence of the helicity on 
the Higgs to W mass ratio, a feature also observed in the generation of Chern-Simons 
number [75, 74]. In the remaining of this section we will focus on results for this 
particular value of the mass ratio. Comments upon the dependence on m^jm^ are 
deferred to section 6.2. 

The late time evolution of the integrated magnetic, longitudinal and transverse 
electric energies, for = 3m^ is presented in Fig. 10. A large fraction of the 
electromagnetic fields generated after SSB is preserved by the time evolution. From 
mt ~ 60 onwards, the transverse energy densities increase with time, again with a 
power law dependence: t", with a = 0.350(1) and 0.330(1) for electric and magnetic 
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Figure 10: We display the log-log plot of the time evolution of the electric (transverse 
and longitudinal) and magnetic energy densities. The exponents of the power law fits 
are: Transverse electric field: 0.350(1); Magnetic field: 0.330(1) and Longitudinal electric 
field: —0.234(2). For = 3m^, ma = 0.42 and p^^^ = 0.15 m, averaged over 200 
configurations. 



energy densities respectively. At these late times, transverse electromagnetic fields 
are composed of an admixture of radiation and long range seed fields. In section 5.2 
we will see how to separate these two components by analyzing the electromagnetic 
field power spectra. Note also that there is a significant fraction of longitudinal 
electric fields, even at the later stages of the evolution. As already mentioned, the 
slow screening of the longitudinal component of the electric field is tied to the presence 
of large charged lumps around magnetic field lines, see Figs. 7 and 8, which persist 
even at late times. 

5.2 Electromagnetic field spectrum 

To investigate whether inverse cascade is active during the late time evolution we 
have analyzed the electromagnetic Fourier spectrum. Figs. 11 and 12 display the 
time evolution of (P|E(/c)|2)/V and {k^\B{k)\^) /V, where E{k) and B{k) are the 
Fourier components of the electromagnetic fields and V is the physical volume. The 
most remarkable feature in the spectrum is the peak at small momenta that develops 
with time, which is distinctly separated from the high momentum component. This 
behaviour suggests that the spectrum contains two uncorrelated distributions which 
describe respectively electromagnetic radiation and the long range electric and mag- 
netic seed fields. Following this indication, we have performed fits to the spectrum 
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Figure 11: We plot {k'^\E{k)'\^) /V vs k, averaged over 150 configurations. The lines 
represent fits to the radiation and seed field electromagnetic components according to Eqs. 
(5.3), (5.5) respectively. Results are presented at mt = 105, 145, 185 and 265. In all cases 
rrijj = 3m^, ma = 0.42 and Pn,;,, = 0.15 m. 



where this separation is made exphcit: 

F{k) = F^""'^(A;) + F'^^^ik) (5.1) 

with F = E OT B. For the expectation values of the electric and magnetic correlators 
we obtain accordingly: 

{\E{k)\^) = ilE^'^'^'^ik)]^) + ilE'^'^m'') (5.2) 

{\m\') = {\B-''im + {\B-\k)\') 

In the remaining of this section we will describe these two components, starting 
with the electromagnetic radiation and ending with the infrared component which 
describes the magnetic field seed. 

5.2.1 Electromagnetic radiation 

The radiation component dominates the electromagnetic energy density, its contri- 




Figure 12: The same as in Fig. 11 but for the magnetic component: (/c^|i?(A:)p)/V. 



bution being a factor of 5-10 larger than the one coming from seed fields. Its profile 
is very well described by: 

= -^r^^, (5.3) 



with we{b) = y + i^E(B) parameters given in Table 2. As illustrated in figures 
11 and 12, this distribution fits very well the high momentum part of the spectrum 
but fails in reproducing the low momentum peak. Eq. (5.3) represents free massive 
thermal radiation with non zero chemical potential at temperatures slightly rising 
with time, which we interpret as an effect induced by the plasma of the ly-fields. 

Similar information can be extracted from the distribution of local values of the 
norm of the transverse electric and magnetic fields. For free photons this should 
follow a Maxwellian distribution (see Appendix C): 

[2/ 3 \3/2 3b2 
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mt 


TE/m 


mE/m 


fiE/m 


TB/m 


mB/m 


fiB/m 


105 


0.32(1) 


0.77(1) 


0.61(1) 


0.32(1) 


0.66(1) 


0.60(1) 


125 


0.33(1) 


0.74(1) 


0.58(1) 


0.33(1) 


0.61(2) 


0.57(2) 


145 


0.34(1) 


0.75(1) 


0.58(1) 


0.33(1) 


0.60(2) 


0.56(2) 


165 


0.34(1) 


0.76(2) 


0.59(1) 


0.34(1) 


0.61(2) 


0.57(2) 


185 


0.34(1) 


0.82(1) 


0.63(1) 


0.34(1) 


0.65(2) 


0.60(2) 


205 


0.35(1) 


0.84(1) 


0.64(1) 


0.34(1) 


0.64(2) 


0.59(2) 


245 


0.35(1) 


0.93(1) 


0.68(1) 


0.35(1) 


0.64(1) 


0.59(2) 


265 


0.36(1) 


0.93(1) 


0.67(1) 


0.35(1) 


0.65(2) 


0.59(2) 



Table 2: Parameters of the fit to the high momentum part of the transverse electric and 
magnetic spectra in Eq. (5.3), for = 3m^, ma = 0.42 and p^,^ = 0.15m. Errors in 
parenthesis combine both systematic and statistical effects. 



where B = \B{x)\. Our data does indeed reproduce this behaviour at late times. 
In Fig. 13 we display the time evolution of the distribution of magnetic field norms, 
starting from mt = 5. Although initially the distribution differs significantly from the 
Maxwellian one, it is approached as time evolves and photons thermalise. There is, 
however, a systematic mismatch when we fit the tail of the Maxwellian distribution, 
even at large values of mt. This signals again a deviation from free radiation, like 
the one observed in the low momentum part of the magnetic and electric spectra. It 
is in this deviation where the contribution of the seed magnetic fields reside. 
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Figure 13: We show the time evolution of the distribution of magnetic field norms. Left: 
For = 3m^ we display the log oi P[B) / B"^ vs B"^ jB'^^-^ (i.e. normalized to the value at 
the peak of the distribution) . Right: For = 3m^ we compare the initial distribution 
of the local magnitude of the magnetic field at mt = 5 with the one obtained at mt = 265, 
the latter fitted to a Maxwellian distribution. The fit to the mt = 5 data is described in 
Appendix D. 
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5.2.2 Electric and magnetic seeds 

We turn now to the analysis of the infrared part of the spectrum, which is the relevant 
one for the generation of the LSMF seed field. This low momentum part has been 
fitted to: 

1 ^ 9k- 

i-ilE'^^'^k)]') = . \ , (5.5) 



V Q0B{wB-ij.B) — \ 



with we{b) = Y ~ + ^1(5) parameters given in Tables 3, 4. This 

could represent again massive radiation at non-zero chemical potential if it were not 
for the peculiar shift /cq in the frequency w. We interpret the value of ~ 0.3m as 
a characteristic momentum scale of the long range electromagnetic fields. 



mt 


TE/m 


TTlE/m 


fiE/m 


k%/m 


105 


0.11(1) 


0.33(5) 


0.30 (4) 


0.29(1) 


125 


0.13(1) 


0.24(4) 


0.22(3) 


0.29(1) 


145 


0.14(1) 


0.21(5) 


0.18(3) 


0.30(1) 


165 


0.13(1) 


0.25(5) 


0.23(3) 


0.29(1) 


185 


0.09(2) 


0.49(8) 


0.48(6) 


0.27(1) 


205 


0.11(1) 


0.36(6) 


0.35(3) 


0.29(1) 


225 


0.10(2) 


0.39(10) 


0.38(3) 


0.28(1) 


245 


0.11(1) 


0.37(7) 


0.35(3) 


0.30(1) 


265 


0.10(1) 


0.45(7) 


0.44(4) 


0.28(1) 



Table 3: Parameters of the fit to the low momentum part of the transverse electric 
spectrum in Eq. (5.5), for = 3m^, ma = 0.42 and p^^^-, = 0.15 m. 



A quantitative estimate of the energy density and correlation length of the seed 
electromagnetic fields can be obtained from our fits to the low momentum part of 
the spectrum. The mean energy density is computed from the integral of the seed 
field spectrum as 

k 

with F = E{B). The correlation length, ^e{b), is extracted from 

^ = _ with k^ = ^ i ^—^ . 5.7 

k Ek\F-m\' 

Table 5 and Fig. 14 summarise our results. We have tested finite volume indepen- 
dence by comparing two different physical volumes: = 0. 125 m and p^-^i^ = 0.15 m. 
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mt 


Ts/m 


iriB/m 


jiB/m 


hPg/m 


105 


0.11(1) 


0.32(7) 


0.30(3) 


0.29(1) 


125 


0.13(1) 


0.24(7) 


0.21(4) 


0.31(1) 


145 


0.13(1) 


0.24(6) 


0.22(3) 


0.29(1) 


165 


0.13(1) 


0.27(6) 


0.23(4) 


0.29(1) 


185 


0.13(2) 


0.18(10) 


0.16(8) 


0.32(3) 


205 


0.11(1) 


0.31(7) 


0.29(4) 


0.30(1) 


225 


0.11(1) 


0.26(5) 


0.25(4) 


0.31(1) 


245 


0.10(1) 


0.37(9) 


0.36(2) 


0.29(1) 


265 


0.11(2) 


0.33(9) 


0.32(3) 


0.30(1) 



Table 4: Parameters of the fit to the low momentum part of the magnetic spectrum in 
Eq. (5.5). For = 3m^, ma = 0.42 and p^in = 0.15 m. 

The numbers in Table 5 come from an average of the results obtained at these two 
physical volumes, with errors given by the dispersion between them. 

We obtain a magnetic seed whose mean energy density increases linearly with 
time. Within the time ranges we have analysed, its fraction to the total comes out 
to be of order ~ 10~^. Assuming the magnetic field expands as radiation, this would 
give magnetic fields today of order 0.5 /iG, which are in the range of the observed ones 
in galaxies, and even in clusters of galaxies, where no-extra amplification through a 
dynamo mechanisms is expected. 



mt 


(Psid)(xl02) 


m^E 


(Psld)(xl02) 


m^B 


105 


0.62(5) 


25.3(1) 


0.58(3) 


25.7(6) 


125 


0.73(2) 


25.2(1) 


0.61(1) 


24.5(9) 


145 


0.76(4) 


24.8(9) 


0.72(2) 


24.8(3) 


165 


0.76(4) 


26.0(10) 


0.77(1) 


25.4(6) 


185 


0.83(1) 


27.6(1) 


0.79(2) 


26.0(10) 


205 


0.89(2) 


27.7(2) 


0.79(6) 


27.2(5) 


225 


0.91(5) 


27.9(5) 


0.87(1) 


28.0(5) 


245 


1.06(9) 


27.6(4) 


0.88(1) 


28.1(2) 


265 


1.12(7) 


27.9(2) 


0.92(2) 


28.4(7) 



Table 5: Fraction of total energy and correlation length of the seed electromagnetic fields. 
They are both derived from the infrared spectrum as described in Eqs. (5.6) and (5.7). 
The results are obtained by averaging (over 150 configurations) the values obtained for 
Pmin = 0.15 m and Pmin = 0.125 m, with errors reflecting the dispersion between them. 
Data correspond to m^ = 3m^, ma = 0.42. 
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Figure 14: We show the time evolution of p^^^^ (left) and m^B (right), for = 
3m^, ma = 0.42. The results are obtained by averaging the values obtained for 
PmiiL = 0.15 m and Pmjn = 0.125 m, with bands representing the dispersion in the errors. 
The fits are /J^ed/Po = 0.0035(5) + 2.3(3) x lO'^mt and m^s = 20.1(4) + 0.033(2)mf 
respectively. 

Concerning the correlation length, it is difficult to make a definitive statement 
about the presence of inverse cascade, given the small time scales we can explore with 
our numerical simulation. Nevertheless, within the time span we have analyzed, our 
results clearly show a linear increase of the magnetic correlation length with time 
(see Fig. 14). This result is robust under changes of p^^,-, and lattice spacing. The 
observed growth is described by m^Bif) = 20.1(4) + 0.033(2)mt, giving at mt = 265 
a characteristic length scale for seed magnetic fields of order m^Bi'fn't = 265) ~ 
30(1). This is much larger than the thermal correlation length, m^thermai ~ 10, and 
represents a significant fraction of the physical volume. It also implies a considerable 
increase from the initial value at mt = 5, obtained from the initial spectrum to be 
m^Bi^nt = 5) ~ 17. From these results we can safely conclude that the time evolution 
has succeeded in amplifying the correlation length of the magnetic seed generated at 
SSB. Nevertheless, a more detailed study, including plasma effects, would be required 
to determine whether ^ will be further amplified at late times. 

In addition to the direct analysis of the spectrum we have also followed an al- 
ternative strategy to separate both the magnitude and the scale of the magnetic 
remnant from the radiation bath. A common way to do this, which has been exten- 
sively used in the literature, is through the computation of several spatial averages of 
the electromagnetic fields. Following Ref. [28], we introduce the following averages: 

• A line average: 





with C a straight line of length /. 
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Figure 15: We show the dependence with ml of the three spatial averages (5.8) — (5.10), 
for mt = 245. The hnes are extracted from our fits to the infrared and radiation parts of 
the spectrum. Note that the fall-off at large distances is just a volume effect. 



The average magnetic flux over a surface of area P: 

1 



i?(2)(0 = 7^ B-dS, (5.9) 



A volume average: 



5(3)(/) = ^ / Mx. (5.10) 



P 



s 



As discussed in Ref. [28], the spatial and statistical averages {B'^^{1)) can be easily 
computed in terms of the spectra of the fields. For instance, the line average for a 
volume V is given by: 

(^(i)(0) = yE^^'(^i'0 (5-11) 

k 

with 

W(k.,,) (5.12) 

Analogous expressions can be found for the other two quantities. The advantage of 
these averages is that they filter out the high momentum part of the spectrum and 
allow to recover, at large /, information about the low momentum modes. We have 
checked that our fits to the spectrum correctly reproduce the spatial averages. This 
is illustrated in Fig. 15, where we present results for the three averages at mt = 245 
compared with the predictions obtained from our fits to the spectrum. The quality 
of the agreement can be considered very good given that the continuum lines are 
directly obtained from the fits to the spectrum (Eqs. (5.3), (5.5) and Tables 2-4), 
and not result of a fit to the spatial averages. 

To summarize, we have found evidence of the presence of a long range helical 
magnetic field, whose amplitude and correlation length are linearly increasing with 
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Figure 16: Comparison of the fraction of total energy carried by electric (transverse and 
longitudinal) and magnetic fields. Top left: for two different values of the minimum mo- 
mentum: Pnjin = 0.1 and 0.15 for fixed ma = 0.65. Top Right and down: 3 different lattice 
spacings ma= 0.65, 0.52,0.42, for the longitudinal, transverse and magnetic components of 
the energy. The lines are the extrapolation of the results to the continuum a ^ limit. For 
rrijj = 2m^ which, from the point of view of lattice artefacts, is the worst case situation. 



time. This is accompanied by the growth of a similar long range electric field. The 
fate of these electromagnetic field depends on the subsequent evolution of the plasma 
which is not addressable within our classical approximation and would require a 
magnetohydrodynamics treatment including the effects of fermion fields. Our results 
for the power spectrum of the seed fields can be used as initial conditions for a MHD 
treatment as the one developed in Ref. [31]. 

6. Dependence on methodological and model parameters 

In this section we study the (in-) sensitivity of our results to the lattice and finite 
volume artefacts. We conclude that all our qualitative results are unaffected by both 
types of approximations. Furthermore, we estimate the size of the systematic errors 
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Figure 17: Left: Lattice spacing dependence of the magnetic susceptibility for mn = 
2m]y, ma = 0.65, 0, 52, 0.42 and = 64, 80, 100. Right: Temporal lattice spacing depen- 
dence of the magnetic susceptibility for mat = 0.05,0.025. 



induced by these cut-offs. The lattice artefacts, though sizable, follow the expected 
O(a^) dependence allowing an extrapolation of the most relevant quantities to the 
continuum limit. 

We also analyze the dependence of our magnetic field production mechanism 
on the Higgs to ly-boson mass ratio m^/m^. It follows from our scenario that, 
initially, the helical susceptibihty xh is independent of the Higgs self-coupling. At 
later times however, we observe a non-monotonic dependence upon the mass ratio, 
which is maximal at our intermediate value m„/m,„ = 3. 

H / W 

6.1 Lattice and finite volume artefacts 

In order to determine the size of the errors introduced by our numerical approach, 
we have performed simulations at different values of the physical volume and of the 
spatial and temporal lattice spacings. The list of simulation parameters is given 
in Table 1. The selection of values implies a delicate compromise among different 
factors. As shown in Ref. [73], to avoid important finite volume effects, we need 
lattices with momentum discretization p^^^ = 27i/L < 0.15 m.. On the other hand, 
concerning lattice artefacts, we have seen in Ref. [74] that cut-off independence of 
certain particular quantities (as the Chern-Simons number) requires m^a < 0.3. 
Most of our lattices satisfy both requirements. 

In Figs. 16 and 17 we present results exhibiting the lattice and finite volume de- 
pendence of the electromagnetic energy densities and of the magnetic helicity. They 
correspond to the most disfavourable case of = 2m^. No noticeable dependence 
on the volume is appreciated. Lattice spacing artefacts are somewhat stronger but 
do not change the general pattern of behaviour. To analyse this effect in more detail, 
we display in Fig. 18 the dependence of the electromagnetic field energy densities 
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Figure 18: Continuum extrapolation of the magnetic, transverse electric, longitudinal 
electric and Z-boson susceptibility. For = 2m^ and mt = 95, 145, 190, and mt = 
60, 190 for the susceptibility. 



and Z-boson susceptibility at various times. In all cases the results are consistent 
with the expected quadratic dependence. This allows the extrapolation of the results 
to the continuum limit, displayed as a continuous line in Figs. 16 and 17. The right- 
hand side of the last figure shows that for the case of the magnetic susceptibility 
the values obtained for the different lattice spacings are compatible within statistical 
errors. Nonetheless, assuming that the lattice spacing dependence depends smoothly 
on time, we can obtain an extrapolation to the continuum limit lying approximately 
5% above the values obtained for the smaller spacing. 

With respect to finite size effects, long range quantities are expected to be the 
most affected. Thus, it is essential to test that the low momentum part of the 
magnetic power spectrum is not biased by finite volume artefacts. In Fig. 19 we 
present results for p^^^ = 0.125 m and 0.15 m. The agreement is very good for the 
ratio m^/m^ = 3 and preserves the same quality for the other 2 values of the mjj to 
m,,, mass ratios that we have studied. 
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Figure 19: We plot {k'^\B{k)\'^) /V vs k for the magnetic component of tlie electromagnetic 
energy. A comparison is made between results at p^in = 0.125 m and p^in = 0.15 m. Results 
are presented at mt = 105 (Left) and 265 (Right). For = 3m^ and ma = 0.52. 



6.2 The Higgs to W boson mass ratio 

Most of the resuhs presented in the previous sections correspond to a Higgs to W- 
mass ratio of 3. Quahtatively the picture remains the same for the other two ratios 
analyzed: = 2 and = 4.65 m^. In Fig. 20 we compare the electromagnetic 
energy densities and helical susceptibility as a function of time for different values 
of the ratio m^/m^. We have chosen here not to normalize the energy densities to 
the total one, in order to exhibit the independence of the initial magnitude of the 
electromagnetic fields and helical susceptibility on the value of Higgs self-coupling A, 
which also determines the mass ratio. Other features of the initial configuration such 
as string lengths and widths are also A- independent, and depend only on the mass 
parameter M that fixes the Higgs Gaussian random field (see Appendix D). This 
A-independence is preserved in the first Higgs oscillation but lost afterwards, once 
non-linearities and the presence of the gauge fields modifies the dynamics. At late 
times equipartition would indicate that the total fraction of energy density carried 
by the electromagnetic field would again become A-in dependent. Since po = 'm^/4:X, 
the fraction of energy densities in units of should tend to behave as 1/A at late 
times. This is indeed the tendency observed in the data. 

7. Conclusions 

In this paper we have analyzed the production of primordial magnetic fields in a 
model of low-scale EW hybrid inflation. Some partial aspects of our study were 
anticipated in Ref. [78]. For that purpose we have studied, with the help of lattice 
non-perturbative techniques, the preheating and early reheating periods after the 
end of a inflationary period. Our work includes, for the first time, the full Standard 
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mt mt 

Figure 20: Time evolution of the energy densities, in m'^ units, in the: Top Left: lon- 
gitudinal electric field; Top Right: transverse electric field; Bottom Left: magnetic field. 
Bottom Right: XH in "i^ units. Energy densities are not normalized to the total energy 
density in order to emphasize A independence in the initial stages of the evolution. 



Model, SU(2) U(l), gauge degrees of freedom. The period of low-scale inflation 
which sets the initial conditions of our work could be brief. We do not need the full 
60 e-folds that are necessary to account for the CMB anisotropics. All that is needed 
is a period of thermal inflation at the EW scale which would cool down the universe 
during at least 10 e-folds, and set the stage for a cold (quantum) EW transition. 
The metric fluctuations responsible for large scale structure could be produced at 
the primordial (high energy scale) inflation. This secondary stage only redshifts scales 
by another factor, but is irrelevant for horizon size fluctuations today, while is 
enough to erase all relativistic and non-relativistic species. This scenario was flrst 
proposed in Ref. [70] and has recently been considered in Ref. [69]. 

The main results of our work can be summarized in the following three observa- 
tions. First, this set up provides a concrete realization of the mechanism proposed 
by Vachaspati [34] and Cornwall [38], by which inhomogeneities of the Higgs fleld 
phases act as sources for the generation of magnetic flelds and — this is essential 
— with non-trivial helicity. To the best of our knowledge this is the flrst time that 
this mechanism has been observed in a fully non-perturbative set-up. Second, the 
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generated magnetic field would have, when red-shifted until today, an amplitude of 
~ 0.5 fiG. This is enough to explain the values of magnetic fields observed in clusters, 
while those in galaxies would require a small amount of enhancement via the usual 
dynamo mechanism. Third, the correlation length of the generated magnetic field 
grows linearly with time within the time span we have analyzed. For = 3m^ we 
find mC,B ~ 0.03 mt, as shown in Fig. 14. This linear growth seems to be sustained 
by the non-trivial dynamics of the plasma made of Pt^-bosons and could be expected 
to hold until the decay of the Higgs, the W and the Z bosons into light fermions. 
Our approach does not allow us to extrapolate these results from then onwards. 
Nevertheless, the helical nature of the generated magnetic field warrants that the 
effect of the primordial plasma would be that of preserving and even amplifying the 
magnitude of the helicity and the magnetic field correlation length [31]- [33]. 

We have distinguished three different stages in the evolution after inflation ends: 
tachyonic growth of the Higgs-fleld low momentum modes, symmetry breaking and 
late time evolution after SSB. In what follows we will summarize the main features 
characterizing each of these stages. 

During the flrst tachyonic stage, non-linearities in the Higgs potential and gauge 
fields can be neglected and the quantum evolution of the system can be exactly 
solved. Quantum fiuctuations of the Higgs-field infrared modes are described by a 
multi-component Gaussian random field. As described in detail in section 4, mag- 
netic fields are already present at this stage with a non-trivial helical susceptibility 
directly related to the winding number susceptibility of the Higgs as a Gaussian ran- 
dom field. Although SU(2)® U(l) gauge fields are very small at the end of infiation, 
the magnetic fields arise through the presence of inhomogeneities in the Higgs field 
phase, thus corroborating Vachaspati's conjecture. Along this period, the spatial dis- 
tribution of the magnetic field is determined by that of the Higgs field, a feature that 
is maintained and even enhanced during the second stage of evolution corresponding 
to symmetry breaking. 

The period of SSB arises via the formation of bubbles in the Higgs field norm 
that expand with time and collide with each other. Magnetic fields are squeezed 
by the expansion in string-like structures localized in the regions between bubbles 
(see Fig. 4). This stringy structure is reproduced both in the helicity density and 
in the Z boson magnetic field density. We have estimated a characteristic string 
separation during this period of m^l ~ 14. Linked to the appearance of the mag- 
netic strings we find a non trivial distribution of electric fields and ly-boson charge 
and current densities. Most remarkably, we see a very non-trivial distribution of 
the charge density with the formation of extended charged clusters which track the 
position of the magnetic string. This separation of unequal charges induces electric 
fields in the plasma. We observe both transverse and longitudinal electric fields also 
correlated with the string locations. The clusters persist for a very long time and, 
as a consequence, we observe a very slow screening of the longitudinal electric field 
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with time, see Fig. 10. We conjecture that these electric fields will be erased as soon 
as the plasma of IV-bosons decays into light fermions moving close to the speed of 
light, which will neutralize much faster than the heavy H^-charges. 

The third stage of evolution after SSB is characterized by a very slow approach 
to thermalisation. To claim a feasible mechanism for magnetogenesis we have to 
guarantee that the initial helical magnetic seed is not removed with time. We 
have shown in section 5 that the magnitude of the helical susceptibility grows with 
time with a power-law behaviour, xh oc t", with a = 0.7(1), 0.8(1) and 0.3(1) for 
m^/m^ = 2, 3, 4.65 respectively. At the same time the Z-boson helical susceptibil- 
ity decays also with a power law dependence with time. We have observed that the 
magnitude of the generated magnetic susceptibility does not depend monotonically 
on the Higgs- to W^-mass ratio. Of the values we have analyzed, m^/m^ = 3 is the 
one that generates larger helical fields. 

In order to extract the late time behaviour of the amplitude and correlation 
length of the magnetic field seed, we have performed a detailed analysis of the mag- 
netic field Fourier spectrum for m^/m^ = 3. It shows two well differentiated and 
uncorrelated components: an ultraviolet radiation sector and an infrared peak whose 
amplitude increases with time (see Fig. 12). The radiation tail is well described by a 
Bose-Einstein distribution of massive photons with a non-trivial chemical potential 
at temperatures T ~ 0.23 slowly rising with time. The low momentum part of 
the spectrum carries a fraction / ~ 10~^ of the total energy density. As mentioned 
before, both its amplitude as its correlation length are linearly growing with time 
within the analyzed time span, showing indications of an inverse cascade towards 
the infrared. However, our time scales are not long enough to demonstrate that 
inverse cascade will be sustained at even later times when the composition of the 
plasma changes significantly. For the moment we can, nevertheless, rely on the re- 
sults in Refs. [31]- [33] which show that helical fields are optimally amplified by MHD 
evolution. 

In summary, hybrid preheating at the EW scale could be responsible for the ob- 
served magnetic fields associated with large scale structures like galaxies and clusters 
of galaxies. Both the magnitude and correlation length could be derived from the 
highly non-linear and non-perturbative evolution after EW symmetry breaking. Our 
analysis provides a concrete realization of the mechanism proposed by Vachaspati 
and Cornwall many years ago. This primordial plasma enters a regime in which 
helical magnetic field lines experience an inverse cascade towards larger scales. We 
observe how both their energy density and correlation length grow linearly with time. 
Showing that these magnetic fields evolve as described in the introduction until pho- 
ton decoupling would require a detailed follow up with MHD simulations with initial 
conditions provided by our work. This result would support our proposal that the 
helical magnetic fields produced at the cold EW transition are responsible for the 
observed magnetic fields in galaxies and clusters of galaxies. 
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A. The Lattice Equations of Motion 

To solve the classical equations of motion we discretise them on a lattice preserving 
full gauge invariance. In this appendix we introduce the lattice notation and derive 
the lattice equations of motion for our particular problem. 

As usual, the lattice points are labeled by a vector of integers n — {no,n) in 
terms of which the space-time positions are given by: x — {no at, n a), with and 
a the temporal and spatial lattice spacings related by k = at /a. The adimensional 
lattice scalar fields are derived from the continuum ones as: $l(^^) = a^{x/a) and 
Xiin) = ax{x/a). In what follows we will omit the subscript L, since all fields 
will be lattice fields unless explicitly indicated. The Higgs field is expanded as: 
$(n) = ^^0"(n)cTQ, in the basis of 2 x 2 SU(2) matrices: a = (1, if), with Ta 
the Pauh matrices and 0"^ real coefficients. The Standard Model Higgs doublet is 
obtained trough the projection: (/? = $ (1, 0)^. Gauge fields are given in terms of hnk 
variables: U^{n) and B^[n) for SU(2) and hypercharge fields respectively. They are 
both 2x2 SU(2) matrices, with the peculiarity that B^{n) is diagonal. Expanded 
in the a basis, they read: 

Ki^)^a, BM=Y. W^-, (A.i) 

a=0,--,3 a=0,3 

with u'^ and &^ real coefficients. The continuum limit of the gauge links is as usual: 

f/^(n) ~ et w.4«r„ ^ (A.2) 

where there is no implicit sum in the /i index and where the vector = {at, a, a, a}. 



-39- 



With the previous conventions, the usual U(l) hypercharge transformation is 
implemented by acting on the $ field with right multiplication by a diagonal SU(2) 
matrix: 

cp'{n) = e^"(") ip{n) = $(n) e^'^^^^^^ . (A.3) 

The complete SU(2)® U(l) gauge transformation for the Higgs field then reads: 

$(n) ^ Q(n)$(n)A(n) , (A.4) 

where A(n) = exp{ia{n)T3) represents the U(l) gauge transformation and fl{n) — 
o-Q the SU(2) one. The corresponding transformations of the gauge links 

are: 

Uf,{n) n{n) U^in) + /x) , (A.5) 

where /■}, is the unit vector in the n direction. 

It is useful to introduce a lattice covariant derivative operator defined by: 

{D^^)in) = U^in) $(n + fi) B^{n) - $(n) , (A.6) 

and its adjoint: 

{D^^){n) = C/t(n - jl) - jl) Bl{n - jl) - $(n) . (A.7) 

In addition we introduce forward and backward ordinary lattice derivatives given by: 

(A^/)(n) = /(n + /i)-/(n), (A.8) 
(A^/)(n) = /(n-/i)-/(n). (A.9) 

The discretization of the pure gauge part of the Lagrangian is done in terms of 
the plaquette fields: 

PA'n) = U^in) U,in + /i) U\n + z>) U\n) , (A. 10) 

P;t{n) = B^{n) B,{n + fi) B\n + D) B\n) , 

with the transformation properties: 

P^,(n) ^n{n)P^,{n)n\n), (A.ll) 
P^,tin) - A(n) P,t(n) A\n) = P;l{n) . 

The pure gauge discretized Lagrangian then reads: 

i 9y 
W)H = -\Y.^A^-Po^{n)]-^Y.^I[l-P^,{n)]. (A.13) 
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And the complete lattice Lagrangian is: 



LL(n) = LyH + Lsu(2)(n) + Tr|(D^$)t(n) (D^<l>)(n)} (A.14) 



+ ^A^xH A^x(^) - V{^{n), x{n)) 



where all the derivatives are lattice derivatives and all matter fields are adimensional 
lattice fields. To simplify notation we have introduced the lattice metric tensor rj^^ 
with non-zero elements: rf^ = I/k and 17" = — /t, i = 1,2,3. This allows to raise 
four-dimensional indices in the usual way. The potential, l^($(n), has the 

explicit form: 

xH) = - Ml Tr{$t(n)$(n)} + A (Tr{<l>"f(n)$(n)})^ (A.15) 



+ ^x'H+^7VHTr{<l>t(n)<|.(n)}, 

with Ml = am, fiL = a ^, and where g and A are the same coupling constants 
appearing in the continuum Lagrangian. 

We have now all the necessary ingredients to write the lattice equations of motion. 
They are derived by imposing that the variation of the lattice action with respect to 
each of the fields in the Lagrangian vanishes. We obtain: 



(A.16) 



(A^A^X)H = «:|/ii + 2(7'Tr[<l>t(n)<l>(n)]|x(n), 
{D^D''^){n) = /t| -M^ + /x^H + 2ATr[$t(n)$(n)]}$(n), 

{Dl Fn{n) = K Ji^{n) , 



with the currents given by: 



2 
2 



- $t(n) (L"^$)(n) ^3/2 
J 3 



(A.17) 



where the sub-index 3 in the second equation denotes the component, of the term 
between brackets, along CT3/2. The covariant derivatives, (D^$)(n) and (i)^$)(n), 
are given by Eqs. (A. 6), (A. 7). We have also introduced two additional covariant 
derivative operators: and D^, obtained from the standard one by setting either 
the hypercharge or the SU(2) gauge links to the identity, i.e.: 



(D;^$)(r2) = U^{n) $(n + fi) - $(n) 



(A.18) 
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The corresponding expressions for the plaquette fields are: 

(Z}^Pp,)(n) = U,{n) Pp,{n + /i) f/t(n) - P,,(n) , (A.19) 

The tensors G^j, and F^^,, appearing in the equations of motion, are defined from the 
traceless part of the plaquettes by: 

P^in)] , (A.20) 
PuM] ■ 

In order to simphfy the problem of solving the lattice equations of motion it is 
convenient to fix the temporal gauge, realized on the lattice by fixing the temporal 
component of the hypercharge and SU(2) links to unity: BqIu) = 11 , Uo{n) = 11. In 
this gauge, the lattice equations of motion can be used to solve for the fields at time 
^0 + 2 in terms of the fields at times uq and uq + I. The lattice equations associated 
to the gauge fixed degrees of freedom become constraint equations analogous to the 
continuum Gauss law: 

{DtG'')in) = Kj'in), (A.21) 

As proved in Ref. [74], these constraints are preserved by the lattice evolution. It 
is hence sufficient to impose them on the initial conditions. The way this is done 
for our numerical simulations follows exactly the procedure described in Ref. [74] for 
SU(2) where we refer the reader for further details. 

B. Lattice version of the Maxwell equations 

In this appendix we present the derivation of the lattice version of the Maxwell 
equations used in order to define the W charge and current densities. Starting from 
the continuum expressions: 

VE{x) = p{x) , Vfix) + dop{x) = (B.l) 
VB(x) = 

V X E{x) + do B{x) = 

V X B{x) -doE{x) =j{x). 

we look for a discretization that preserves the Bianchi identities. 



G 



2gY 

i 



KM 

[P,u{n) 
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In section 3.2, we have defined the electromagnetic lattice field strength, F^j^{n), 
in terms of clover averaged Z and B field strengths. The clover average of a space- 
time tensor, like Foi{n), is given by: 

(Foi(n))ciov = \ {F^^{n) + Fo,(n - 6)) , (B.2) 

while for a spatial tensor we have: 

(i^ijH)ciov = \ (F,j{n) + Fij{n - i) + F,,(n - j) + Fij{n - i - j)) . (B.3) 

From them we extract the lattice electric and magnetic fields: 

Ei{n) = ^— (Fio(n))ciov , Bi{n) = eijk {Fjk{n))ciov ■ (B.4) 

eaat 2ea^ 

^ — * — # 

The electromagnetic E and S fields, defined above, verify the following Bianchi 
identities: 

■ B{n) = , X E{n) + - AoB{n) = , (B.5) 

K 

where we have introduced an improved lattice derivative given by: 

(Ai/)W = ^(/(^ + A)-/(r^-A)). (B.6) 

We now define accordingly the longitudinal and transverse components of the electro- 
magnetic fields. The projection is done in momentum space with Fourier transformed 
fields: 

with lattice momenta: fcj = 27mi/ (Ng a), rii G Transverse components. At, of a 
vector A, are defined such that q ■ At = 0, where: 

q= - (v-v*), with v, = - {e-'^^^ - 1) , (B.8) 

2 Qj 

and with q the unit vector in the direction of q. 

The electromagnetic, Fourier transformed, charge and current densities are com- 
puted through: 

p{k) = q-E{k), (B.9) 

j{k) = q^ B{k)-—KQ E{k). (B.IO) 
at 
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C. Thermal radiation 



In the present appendix, we prove the relation: 

mx)r) = {:\B{x)r:)Q^^T). (C.l) 

where the left side of the equality is calculated using the Maxwellian classical distri- 
bution: 

whereas the right hand side is calculated using the thermal quantum distribution in 
the canonical formalism. Thus, 

MRm|2nA _ Tr(: : p) 

(• \^{X)\ -jQiT) = 7f^^ • (<--^) 

where p is the canonical distribution density matrix: 

p = e-^. (C.4) 



and : O : denotes normal ordering of the operator O. By performing the integral in 
Eq. C.2 we obtain the classical thermal averages: 



(2v -I- IV 

{\B{x)n= ^ >- {{\B{x)\')r. (C.5) 

Our goal is then to compute the thermal quantum averages 

(: {B{x) ■ B{x)T ■■)Qm 

on the canonical distribution at temperature T. The only terms of the normal-ordered 
operator that contribute to the expectation values must be diagonal in momentum 
space. If we single out that part we obtain 

■.{B{x)-B{x)r:=f[(Yl / ^^^X(^'*^)aa.(^l))Ga,...„„(ri,...Q+X, (C.6) 

where X denotes the part that does not contribute to the expectation value and G 
coefficient to be specified later. 
Next we can evaluate the thermal average of the operator part, which can be 
expressed as a product of n{ki,ai), the mean number of photons of momenta ki and 
polarization a^. Hence, we arrive at 

ai...a„ (^1) • • • ■ 



i=l \ a 
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Now we should unfold the form of the coefficient G. It is given by 

where the sum is over all the permutations of the 2n indices, and 

The sum over all permutations follows from taking all creation annihilation operators 
as distinguishable and assigning them to each of the 2n magnetic fields. Nonetheless, 
since we are integrating over all values of momenta one has to divide by n\ to eliminate 
double-counting . 

Now we will introduce the matrix M, given by 

Mij = Y^ / dkn{k, a)vi{k, a)vj{k, a) = Xdij . (C.7) 

The left-hand side is a consequence of rotational invariance. Substituting in the 
previous formulas we get 

~\-^hh ' ' ' ^hn-lhn ^ ] ^'<t{1)'<t(2) ' ' ' ^'<T(2n- 1) '<T (2n) ' 
0'ei>2n 

The sum over permutations can be factored as follows 

cr&S2n pairings pair 

where a pairing is an arrangement of the 2n indices into pairs (equivalently a per- 
mutation made entirely of 2-cyclcs). The rest of the calculation is very much like a 
calculation to nth order in perturbation theory in a model with 2-leg vertices given 
by the M matrix and a propagator given by the identity matrix. All diagrams are 
now characterized by ni, the number of 1-cycles (loops), where 1 runs from 1 to n. 
Applying the standard Feynman rules one arrives at 

The factors 21 and nil provide the order of the symmetry group of the diagram. The 
21 term is associated with cyclic permutations of the vertices and to a change in 
orientation. In the previous formula, the sum over rii runs over all possible integers 
subject to the constraint Irii — n. One can actually perform this sum. Setting 
M — XL our expression becomes proportional to A". Thus, we can eliminate the 
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constraint on the rii by summing over n. The constrained sum can be obtained from 
the unconstrained one by selecting the term proportional to A". Hence, 

where D is the space dimension, which is 3 in our case. This quantity is the generating 
function of all the quantum averages. Differentiating n times with respect to A we 
extract the n-th term that we were looking for: 

(: im ■ B{x)r :)q(t) = (2A)"(2n + 1)!! (C.8) 

The result for D = A = 1, given by serves to crosscheck the result. From the 
previous equation (C.8) we get A = (l/6)(: :)q(t) allowing to re-express eq. 

(C.8) in the form of eq. (C.5). 

To conclude we give the expression of (: |i?(x)p '■)q{t) in terms of the tempera- 
ture. Taking the trace of eq. (C.7) we obtain: 

{■.\B\'{x) :)q(t) = I dkj^^2kn{k,a). 
Taking into account n{k,a) = (e'^/^ — 1)~^, we can perform the integration: 

The sum over n is the known ^(4) = 7r^/90, leading to: 
D. Gaussian Random fields 

In this appendix we revisit the predictions of the Gaussian random field model. As 
explained in Section 3 of the paper, the initial conditions produced by the quantum 
evolution shortly after inflation ends are of this type. Furthermore, this distribution 
seeds the generation of magnetic fields and Chern-Simons number. There is an 
extensive literature (see Refs. [86]- [89]) on Gaussian random fields and some of the 
analytic predictions have been included in our previous papers. However, here we 
are dealing with multicomponent fields and some of the predictions and methodology 
do not hold in this case. Besides, there are many more relevant observables directly 
related to the Physics issues addressed in this paper. To explore these matters within 
this paper, we have felt satisfied with its numerical study. Since gauge fields and 
non-linearities do not play a role at this stage, we have profited to increase statistics 
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and test systematic errors at a low computational cost. These results can then be 
used as a reference to compare with our full-theoretical ones. 

Our Gaussian random field is homogeneous and isotropic. The power spectrum 
was set to match the one produced by the quantum evolution of the Higgs field cou- 
pled to a linearly time-dependent inflaton and neglecting the Higgs self-interaction. 
The details and nomenclature are explained in our previous paper [73]. We recall 
that the Higgs field has 4 real components which are independent random variables 
with identical power spectrum which, for simplicity, is fitted to a simple form which 
reproduces nicely its shape: 

p(k,t) = -Ar^k\A(t)e-''^'^'''^'^' + l)&(VWt^ - k) (D.l) 

where V is the inflaton velocity at the end of inflation, A and B are time-dependent 
parameters and 6 is the Heaviside step function. 

It is interesting to be able to trace the dependence of our results on the different 
parameters that enter our model. Fortunately, this dependence is greatly encoded in 
two scales that characterize the Gaussian random field. One scale fixes the magnitude 
of the Higgs field. We choose this scale to be the dispersion a of the field at one 
spatial point. Notice that the physical scale v, giving the expectation value of the 
Higgs field in the true vacuum, has not yet entered the scene, since the Gaussian 
random field is generated before the self-interaction of the Higgs field affects the 
evolution. It is precisely the comparison between a(t) and v that must be taken 
into account in fixing the range of values of the initial times ti for the subsequent 
non-linear classical evolution of the system. 

In addition, the other scale of the problem is a length scale associated to the 
Gaussian random field as follows: 

1 _ JfP{k,t)e 

e!= jfp{k,t) 

With our choice of velocity V = 0.024 at mti = 5 we obtain a = 0.139 f, for 
TTijj = 3m^. Thus, we are safely in the region where non-linearities are still small. On 
the other hand m^o = 3.09, which determines the adequate ranges of the ultraviolet 
and infrared cut-off of our numerical procedure. At mti = 6.5 these numbers have 
changed to cr = 0.204 and m^o = 2.95 respectively. This observation allows us to 
give results in a way that are valid for all the values of initial times employed in this 
work. 

In line with previous analysis, we will present our results for the density and 
distribution of local maxima of |0|. The density of maxima is given by 0.0140(4) $^q^. 
The distribution of minimum distances among maxima can be studied directly and 
displays an approximate Gaussian distribution with mean 3.1(1) .^o dispersion 
0.62(2)^0- We have also studied the distribution of values of |0| at the maxima, (p. 



(D.2) 



The average height of a peak being 1.52 (4) a. The histogram is much narrower than 
the one obtained for a single component Gaussian random field, and is well-fitted to 
the following expression 



1^°' exp 



2a2 



with a = 10.4(5) and a = 0.44(1) a. Nicely enough the results presented are robust 
as one changes the ultraviolet, infrared cut-offs and time within their safe windows 
(See Fig. 21). Errors quoted are both statistical and systematic. 
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Figure 21: Left: Histogram of peak (local maxima) heights, expressed in cr(tj) units. 
Right: Distribution of the the local magnetic field intensity B = \B{x)\ in (,o{ti) units. 



Now we turn to observables which are characteristic of multi-component Gaus- 
sian random fields. A crucial role is played by the topological susceptibility x which 
is obtained by dividing the mean value of the winding number square by the volume. 
We obtain 1.55(10) x IO'^^q I We can also compute the initial magnetic field dis- 
tribution. Notice that, as explained in the paper, despite the fact that SU(2)xU(l) 
gauge fields are zero at this stage, our formulas induce a non-zero Z field and a 
non-zero magnetic field which is proportional to it. Computing this magnetic field 
at each point of space we obtain a distribution which is well fitted by a formula 

P(S).S'exp{-(|)"}+^e.p{-(|)"j. 

with B = \B{x)\, see Fig. 13. Our best fit values of the parameters are b = 1.89(3), 
ciiCo = 3.0(1) 10-3, /ii = 0.368(3), da^o = 2.61(2), /i2 = 1.34(3) , A = 1.0(5)10"^. 
The initial magnetic field distribution has a slower decrease at large values than the 
Maxwellian distribution obtained at later times. The aforementioned universality 
can be tested here. In particular, it follows that results obtained at different initial 
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times ti should fall in the same curve once normalised by the scales of cr and ^q. This 
is clearly seen in Fig. 21. 

We have also studied the spectrum of the magnetic field to compare it with the 
one obtained once non-linearities set in. In our case the high momentum profile 
differs from the thermal tail displayed at later times. Instead, the high momentum 
tail is well fitted by a function 



where 6^0 = 0.01(1) and c = 0.36(4). 
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